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In this work, we explore asymptotically flat charged thin-shell wormholes of third order Lovelock 
gravity in higher dimensions, taking into account the cut-and-paste technique. Using the generalized 
junction conditions, we determine the energy-momentum tensor of these solutions on the shell, and 
explore the issue of the energy conditions and the amount of normal matter that supports these 
thin-shell wormholes. Our analysis shows that for negative second order and positive third-order 
Lovelock coefficients, there are thin-shell wormhole solutions that respect the weak energy condition. 

In this case, the amount of normal matter increases as the third-order Lovelock coefficient decreases. 

We also find novel solutions which possess specific regions where the energy conditions are satisfied 
for the case of a positive second order and negative third-order Lovelock coefficients. Finally, a 
linear stability analysis in higher dimensions around the static solutions is carried out. Considering 
a specific cold equation of state, we find a wide range of stability regions. 

PACS numbers: 04.20.Jb,04.50.Kd, 04.50.-h 


I. INTRODUCTION 

Wormholes are topologically nontrivial objects which 
connect two separate and distinct spacetime regions 
[1, 2]. A fundamental ingredient in wormhole physics is 
the flaring-out condition of the throat, which in General 
Relativity (GR) entails the violation of the null energy 
condition (NEC). Matter that violates the NEC is de¬ 
noted exotic matter. In fact, wormholes violate all of 
the pointwise and averaged energy condition [2] . Due to 
the problematic nature of these violations, it is useful to 
minimize the exoticity problem of the matter threading 
and sustaining the wormhole. For instance, traversable 
wormholes supported by arbitrarily small quantities of 
exotic matter [3] or by matter respecting the energy con¬ 
ditions [4-7] have been investigated. 

An interesting kind of traversable wormhole is the thin- 
shell variety constructed by surgically grafting together 
two spacetimes to form a geodesically complete manifold 
with a shell placed at the junction interface. In these 
types of wormholes, matter is concentrated on the worm- 
hole throat, which coincides with the junction boundary, 
and thus minimizes the usage of exotic matter. These 
solutions were initially presented in [8, 9] and have been 
extensively investigated in the literature [10, 11]. In GR, 
it was shown that thin-shell wormholes inevitably vio¬ 
late the weak energy condition (WEC) and, consequently, 
various attempts have been undertaken, to alleviate this 
problem, in alternative theories of gravity, for instance, in 
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dilatonic gravity [12], Brans-Dicke theory [13] and DGP 
gravity [14] . An advantage of the thin-shell construction 
is that one may consider a linearized stability analysis 
around the static solution, while preserving the symmet¬ 
ric configuration, and analyse the dynamical evolution 
of the thin-shell [15]. The generalization of the stabil¬ 
ity analysis to the case of higher-dimensional wormholes 
has also been extensively investigated [16]. For instance, 
a linearized stability analysis in the presence of a dila- 
ton, axion, phantom and other types of matter and in 
cylindrical symmetry have also been studied [17]. 

Since string theory proposes higher dimensional grav¬ 
itational objects, it is of interest to explore the possibil¬ 
ity of higher dimensional thin-shell wormhole solutions. 
One of the most interesting higher dimensional alterna¬ 
tive theories of gravity which leads to second order field 
equations for the metric is described by the Lovelock ac¬ 
tion [18]. In this context, thin shell wormholes of second 
order Lovelock theory called Gauss-Bonnet theory in 5- 
dimensions and higher dimensions were extensively anal¬ 
ysed [19]. In the context of third order Lovelock gravity, 
asymptotically flat and charged thin-shell wormholes of 
Lovelock gravity in seven dimensions were constructed, 
using the cut-and-paste technique, and the generalized 
junction conditions were applied in order to calculate the 
energy-momentum tensor of these wormholes on the shell 
[20]. It was found that for negative second order and pos¬ 
itive third order Lovelock coefficients, there are thin-shell 
wormholes that respect the WEC in specific regions. In 
this case, the amount of normal matter decreases as the 
third order Lovelock coefficient increases. For positive 
second and third order Lovelock coefficients, the WEC is 
violated and the amount of exotic matter decreases as the 
charge increases. Furthermore, a linear stability analysis 
was analysed against a symmetry preserving perturba- 
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tion, and the stability of these wormhole were studied. 

In this work, we extend the analysis outlined in [20], 
by providing the explicit form of the surface energy- 
momentum tensor and generalize the analysis to arbi¬ 
trary dimensions. More specifically, we present the ex¬ 
plicit form of the energy-momentum tensor in terms of 
the extrinsic curvature and Riemann tensor for all higher 
dimensions, and find specific solutions, including novel 
solutions which possess specific regions where the energy 
conditions are satisfied for the particular cases of positive 
second order and negative third-order Lovelock coeffi¬ 
cients. Furthermore, we present the generalized junction 
conditions, and analyse the stability analysis for the so¬ 
lutions found, for arbitrary dimensions. Thus, we study 
the effects of the third order Lovelock term in higher di¬ 
mensionality, its role in the satisfaction of the WEC and 
the stability of thin shell wormhole solutions in higher 
dimensions. 

This paper is organised in the following manner: In 
Section II, we present the action of third order Lovelock 
gravity in the presence of an electromagnetic field, and 
construct a thin-shell wormhole in Lovelock gravity, by 
presenting explicitly the surface energy-momentum ten¬ 
sor components on the shell. In particular, to comple¬ 
ment the above Section II, the specific (n+l)-dimensional 
static solutions, that are to be used throughout the work, 
are outlined in Appendix A. In Appendix B, we also in¬ 
troduce the junction conditions in Lovelock gravity, in 
order to construct the thin-shell wormhole. In Section 
III, we analyse the profile of the surface stresses and con¬ 
sider the issue of the energy conditions and the amount 
of normal matter supporting the thin-shell wormhole. In 
Section IV, we study the stability of the thin-shell worm- 
holes under small perturbations preserving the symmetry 
of the wormhole configuration, and explore the dynamical 
evolution of the shell. Finally, in Section V, we conclude. 


II. ACTION AND THIN-SHELL 
CONSTRUCTION 

The action of third order Lovelock gravity in the pres¬ 
ence of an electromagnetic field may be written as 

1 = J d n+1 xy/=g (A + a 2 C 2 + a 3 C 3 - , ( 1 ) 


where a 2 and 03 are the second (Gauss-Bonnet) and 
third order Lovelock coefficients. The term A = R is the 
Einstein-Hilbert Lagrangian, A is the Gauss-Bonnet La- 
grangian, and C 3 is the third order Lovelock Lagrangian, 
which are defined in Appendix A. 

The electromagnetic tensor is defined as F^ u = d^A„ — 
d^A t j. where the vector potential is given by 




(n ~ !) Q rt 
2(n — 2) r n ~ 2 


The (n + l)-dimensional static solution of action (1) is 
given in the form 

dr 2 

ds 2 = - f{r)dt 2 + —— + r 2 (!£];_!, (3) 

/ (r) 

where dH^-i is the metric of a (n — l)-dimensional unit 
sphere and the metric function /(r) is given by 

/(r) = 1 - r 2 ip(r), (4) 

where the solutions for the function ip(r) are outlined in 
the Appendix A. 

To construct a thin-shell wormhole in third order Love¬ 
lock gravity, we use the well-known cut-and-paste tech¬ 
nique. The junction conditions in Lovelock gravity are 
given in Appendix B, where we refer the reader for more 
details. Taking two copies of the asymptotically flat so¬ 
lutions of Lovelock gravity given by Eqs. (3) and (4) and 
solutions of Eq. (A6) and removing from each manifold 
the (n + l)-dimensional region described by 

= {ut < a, a > r h } , (5) 

we are left with two geodesically incomplete manifolds 
with the following timelike hypersurface as boundaries 

T,± = {r± = a,a > r h } ■ ( 6 ) 

Now identifying these two boundaries, £ + = £_ = £, 
we obtain a geodesically complete manifold containing 
the two asymptotically flat regions f2 + and f2_ which are 
connected by a wormhole. The throat of the wormhole 
is located at £ with the following metric 

= — dr 2 + a 2 (r)dfl 2 _ 1 , (7) 


where r is the proper time along the hypersurface £ and 
a(r) is the radius of the throat. In order to study the sta¬ 
bility of the thin-shell wormhole constructions, we allow 
the radius of the throat to be a function of the proper 
time. In addition to this, note that all the matter is 
concentrated on the wormhole throat, £. 

To analyze such a thin-shell configuration, we need to 
use the modified junction conditions. Considering the 
induced coordinates on £, £ a = (r, 0 1 ; i = l...n), the 
extrinsic curvatures associated with the two sides of the 
shell are given by 


where the normal vector n± ( n p n p = 1 ) to the surface £ 
in M is defined by 





d 2 XP 


r p „ 


dXP dX l 


Q£aQ£b ^ Q£a Q£b 


n 


± 

'7 


= ± 



dg og 

dXP dx v 


dg 

dX'T ’ 


( 9 ) 


where g(r,r) is the equation of the boundary £, given 

by 


( 2 ) 


g(r,r ) = r — a(r) = 0 . 


( 10 ) 
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After some algebraic manipulations, in an orthonor¬ 
mal basis {ef,ej; i = 1 ...n — 1 }, the components of the 
extrinsic curvature tensor are given by 


K.% 

r 

A’ 

( 11 ) 

q 

= -*} 
ft J 

( 12 ) 

with 




A = \/a 2 + f (a), 


where the prime and the overdot denote derivatives with 
respect to a and r, respectively. Equations (11) and (12) 
show that the form of the energy-momentum tensor on 


the shell is S a j- = diag (— a,p S l ~), where a is the surface 
energy density and p is the transverse pressure. 

Now using the junction condition (B4), the compo¬ 
nents of the energy-momentum tensor on the shell may 
be written as 

C r _ A / (n- 1 ) 2a 2 (n- 1 ) 

4tt\ a 3 a 3 

x [a’- 3 ( W)] + AAA 

x [15(1 + d 2 ) 2 - 10A 2 (1 + a 2 ) + 3A 4 ] 1, (13) 


P = 


oi 1 f 2r , 2(n — 2)A 
1 87 t \ A a 

4a 2 ( nT , A 3r A , A 3 (n — 4) 


3a 2 

6 A 

a 


{ 3 r A-f (l + a 2 ) + 


(n — 4) 


(l + ft 2 ) + ftC 


2(23 

5a 1 


x J ^ (! + d 2 - A 2 ) 2 + 20dA [3 (1 + d) - A 2 ] 


(n — 6 )A 
5a 


[15(1 + d 2 - A 2 ) 2 


+20 (1 + d 2 ) (3(1+ d 2 )-A 2 )]}} 


(14) 


respectively. From the above equations we see that a and 
p are expressed in terms of the throat radius a(r), the first 
and second derivatives of a(r) and the metric function 
/(a). Note that the surface energy density and transverse 
pressure satisfy the energy conservation equation, given 

by 


iH”- 1 , ) + 4 (“ l ”- 1, )=°- (15) 

The first term in Eq. (15) represents the internal energy 
change of the shell and the second term shows the work 
by internal forces of the shell. 


III. PROFILE OF THE SURFACE 
ENERGY-MOMENTUM ON THE THIN-SHELL 


In this section, we consider the issue of the energy 
conditions and the amount of normal matter that sup¬ 
ports the thin-shell wormhole. The WEC is defined as 
T) lv V^V l/ > 0 for every nonspacelike vector V^. The 
WEC on the thin shell is satisfied provided that a > 0, 
a+p > 0. We restrict our analysis to static configurations 
with a = ao and a = a = 0. For the static configuration, 
Eqs. (13) and (14) reduce to 


CT ° 8 ( 7ra 0 Jo{ 2/o+ 3a 2/o(3 /o) 

+ y~y/o(15 — 10/o + 3/o 2 ) 
oa 0 


(16) 


ao +po 


o-7 t={ ( _2 /o + ao/o) 

87raov Jo *- 

+ -+-7 [3ao/o(l - /o) - 6 / o (3 - / 0 )] 

6d 0 

+ —[3ao/o(l — 2/o + /q) 
a 0 

-2/ o (15-10/o + 3/o 2 )]}, (17) 


respectively, where /o = /(ao) and f' Q = /'(ao). In con¬ 
trast to the case of general relativistic thin-shell worm- 
holes for which ao < 0, thus violating the WEC [16], in 
higher-dimensional thin-shell wormholes one can have, in 
principle, normal matter on the shell. In order to investi¬ 
gate the profile of the surface energy-momentum on the 
thin-shell, we find the amount of matter on the shell, 
which is given by 


F = 


JdrdQ n -i[a 0 S(r - a 0 ) +p r \- 


(18) 


For our case, the shell does not exert radial pressure, i.e., 
p r = 0 , and therefore the amount of matter on the shell 
is 


F = 


2tt2 Oq 1 cto 

r (§) 

4-k% (n — 1 ) 

r(§) 

. 5 3 a[)" _6) 


\/7o 


- a, 


(n—2) 2d 2 aQ 


(n-4) 


(-15+10/o-3/o 2 ) 


-(/o - 3 ) 


(19) 


where fo is the value of /(r) at r = ao- In Einstein grav¬ 
ity, a 2 = d 3 = 0, the matter is exotic both for Reissner- 
Nordstrom (Q ^ 0) and Schwarzschild (Q = 0) thin-shell 
wormholes as one can see from Eq. (19) as ao < 0 (Eq. 
(16)). In Gauss-Bonnet gravity, a 3 = 0 with a 2 < 0, 
one may have thin-shell wormholes supported by normal 
matter [19]. 
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Now, we investigate the condition that thin-shcll 
wormholes may be supported by normal matter in third 
order Lovelock gravity. It is clear from Eq. (19) that the 
sign of F depends on the sign of cr 0 . 


found thin-shell wormholes threaded by normal matter 
on the shell in specific regions. In this case, as depicted 
in Fig. 4, there exist regions for a 0 > 0 (or F > 0) 
and oo + po > 0 and therefore the WEC and NEC are 
satisfied. 


A. Specific case: a 2 = 3d3 = a 2 

For the special case a\ = 3d3 = a 2 , Eq. (16) takes the 
form 

°o = ~ ~ 10q ( q + «g)/o 

+15(a 2 + 2aao + Oq)] ■ (20) 

The sign of the square brackets in the above equation 
determines the sign of a 0 . This is a quadratic equation 
for the metric function fo and one can easily obtain the 
discriminant of the quadratic equation which is equal to 
—80a 2 (a + aQ) 2 . Since the discriminant is negative for all 
values of a, the sign of the brackets is always positive and 
therefore ao < 0. In other words, the matter supporting 
the thin-shell wormhole is exotic for any dimensions. 


B. Specific case: 02 > 0 and as > 0 

For the general solutions of third order Lovelock grav¬ 
ity we first consider the case where cb and a$ are positive. 
In this case, since the metric function is asymptotically 
flat (0 < /o < 1) for a ^ 77 ,, and therefore the factors 
3 — /o and 15 — 10/o + 3/g are positive, the bracket in 
Eq. (16) is always positive. Therefore, the matter on the 
throat is exotic in the case of 012 , £*3 > 0 as in this case 
we have er 0 < 0 (or F < 0 ). 


C. Specific case: 02 < 0 and 03 > 0 

For 03 > 0 and 0:2 < 0, F can be positive and therefore 
the matter is normal. For this case, Figs. 1 and 2 display 
the behaviour of a 0 and eg + Po for n = 8 and 9 respec¬ 
tively. In these plots, there exist regions for cto > 0 (or 
F > 0) and ao+Po > 0 and therefore the WEC and NEC 
are satisfied in this region. Since 0 < /o < 1, the factor 
—15 + 10/o — 3/g in Eq. (19) is negative and therefore 
the amount of normal matter for negative 02 decreases 
as 0:3 > 0 increases, as one can see in Figs. 1-a) and 1-b) 
for ck 3 = 2 and 03 = 1. In Figs. 2-a) and 2-b), we show 
that the amount of normal matter increases as a.i < 0 
decreases. Also, in this case the amount of normal matter 
decreases as the charge q increases (see Fig. 3). 


D. Specific case: ai > 0 and as < 0 

Finally, we can see that for 02 > 0 and a?, < 0, F 
can be positive and therefore, for the first time, we have 


IV. STABILITY ANALYSIS 


In this section, we study the stability of the thin shell 
under small perturbations preserving the symmetry of 
the wormhole configuration. The dynamical evolution 
of the shell results from Eqs. (13) and (14) or by any 
of them and the conservation equation (15). In order 
to analyze the stability and to complete the system we 
use a cold equation of state p = p(a) with r) = dp/da. 
We consider small radial perturbations around a static 
solution with radius ag and obtain the equation of motion 
near the equilibrium solution. In this case, one may write 
p ~ po+rjo(a— (To), where po, ao and rjo are the transverse 
pressure, surface energy density and dp/da at the static 
solution a = a 0 , respectively. Using this linear equation 
of state, then Eq. (15) yields the following solution for 
the surface energy energy 


cr(a) 


/ ap+po \ /ao\ (.n-iy.i+r,) pa 0 _p 0 

V I + 77J V a 7 ' 1 + 77 


( 21 ) 


Substituting the energy density (21) in Eq. (14), leads 
us to the equation of motion, for the radius of the throat, 
given by 


60na 5 a(a) 

(n - 1 ) v / a 2 + /(a) 

3d 3 {8(1 + a 2 ) 2 + 4[1 — /(a)](l + a 2 ) + 3[1 — /(a)] 2 } 
+15a 4 + 10a 2 a 2 [2(1 + d 2 ) + 1 - /(a)] , (22) 


where <r(a) is given in Eq. (21). 

The dynamics of the wormhole throat is determined 
by a single equation with the form of a 2 = —V(a). How¬ 
ever, in third order Lovelock gravity, one encounters a 
fifth order algebraic equation for a 2 and unfortunately 
one cannot obtain an exact result, but, one may perform 
the linearized stability analysis numerically. More specif¬ 
ically, one may solve Eq. (22) for a 2 numerically and ob¬ 
tain the potential V (a) in the equation of motion. Then, 
the wormhole with radius ao is linearly stable provided 
the potential V(a) is minimum at a = ao where a(r) will 
oscillate about ao and stabilize the wormhole. Numerical 
calculations are shown in Figs. 5-7 for n = 8 and n = 9. 
In the analysis, we have chosen the regions ao for which 
matter is normal. As depicted by the figures, for the 
specific parameter ranges, the wormholes are stable for 
negative values for tj and unstable for positive ij = dp/da. 
We emphasize that these results are in agreement with 
those presented in [20], where it was shown through nu¬ 
merical calculations that wormholes are stable provided 
the parameter rj is negative and the throat radius ao is 
adequately chosen. 
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(a) a 3 = 2 


(b) a 3 = 1 


FIG. 1: Specific case: a .2 < 0 and 03 > 0. 00 (continuous) and ao +po (dashed) versus a 0 for a .2 = —2, m = 100, q = 3000 and 
n = 9. These figures depict that for a.2 < 0 and 03 > 0, the WEC can be satisfied. The plots also show that normal matter for 
0:2 < 0 decreases as a 3 > 0 increases. 




(a)c *2 = —2 (b) 02 = —1.5 

FIG. 2: Specific case: 0:2 < 0 and a 3 > 0. ao (continuous) and ao + Po (dashed) versus ao for a 3 = 1, m = 100, q = 3000 and 
n = 8 . These figures depict that for 0:2 < 0 and 0:3 > 0, the WEC can be satisfied. The plots also show that for a 3 > 0, the 
amount of normal matter increases as Q 2 < 0 decreases. 
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FIG. 3: Right: f(r) versus r for ai = —0.6, 03 = 0.2, n = 8 , m = 300 and q = 2000, 3000, 5000 from down to up respectively. 
Left: 10 _4 F versus a 0 for f(r) plotted in the right figure. These figures depict that for 03 > 0 and 02 < 0, the amount of 
normal matter decreases as q increases. 



FIG. 4: Specific case: 0:2 > 0 and 03 < 0. <to (continuous) 
and a 0 + po (dashed) versus a 0 for 02 = 0 . 12 , 03 = — 0 . 1 , 
m = 10, q = 1.4 and n = 8 . 


V. SUMMARY AND DISCUSSION 

Motivated by the fact that string theory proposes 
higher dimensional gravitational objects, we have con¬ 
sidered higher-dimensional thin-shell wormholes in the 
context of third order Lovelock gravity. More specifi¬ 
cally, asymptotically flat charged thin-shell wormholes of 
third order Lovelock gravity in higher dimensions were 


constructed, taking into account the cut-and-paste tech¬ 
nique. We extended the analysis outlined in [20] , by pro¬ 
viding the explicit form of the surface energy-momentum 
tensor in terms of the extrinsic curvature and Riemann 
tensor for all higher dimensions, explored the issue of the 
energy conditions and the amount of normal matter that 
supports these thin-shell wormholes. Our analysis shows 
that for negative second order and positive third-order 
Lovelock coefficients, there are thin-shell wormhole so¬ 
lutions that respect the weak energy condition. In this 
case, the amount of normal matter increases as the third- 
order Lovelock coefficient decreases. We also found novel 
solutions which possess specific regions where the energy 
conditions are satisfied for the case of a positive second 
order and negative third-order Lovelock coefficients. Fur¬ 
thermore, we presented the generalized junction condi¬ 
tions, and analysed the stability analysis for the solutions 
found, in arbitrary dimensions. This enabled us to study 
the effects of the third order Lovelock term in higher di¬ 
mensionality and the satisfaction of the WEC and stabil¬ 
ity of thin shell wormhole solutions. In the latter context, 
a linear stability analysis in higher dimensions around the 
static solutions was carried out, by considering a specific 
cold equation of state, and a wide range of stability re¬ 
gions were found. 
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7.8 7.9 a 8 8.1 7.8 7.9 a 8 8.1 




FIG. 5: V(a) versus a with a 2 = — 2, 0:3 = 2, m = 100, n = 9 and q = 3000 for q = —3 (left) and q = 3 (right). These figures 
depict that wormholes are stable for q = dp/da < 0 and unstable for q > 0 . 

6.6 6.7 6.8 a 6.9 7 6.9 7 7.1 a 7.2 7.3 7.4 
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FIG. 6 : V(a) versus a with 02 = —2, «3 = 1, m = 300, n = 8 and q = 3000 for r? = —2 (left) and 77 = 1 (right). These figures 
depict that wormholes are stable for q = dp/da < 0 and unstable for q > 0 . 



FIG. 7: V(a) versus a with a .2 = 0.12, 03 = —0.1, m = 10, n = 8 and q = 1.4 for rj = —1.5 (left) and q = 1.5 (right). These 
figures depict that wormholes are stable for q = dp/da < 0 and unstable for q > 0 . 









reference IF/00859/2012, and the grants EXPL/FIS- 
AST/1608/2013, PEst-OE/FIS/UI2751/2014 and 

UID/FIS/04434/2013. 


Appendix A: Action and spherically symmetric 
solutions 

We reproduce the action of third order Lovelock grav¬ 
ity, in the presence of an electromagnetic field, here 

1= [ d n+1 x ^/=£ (A + a 2 C 2 + a 3 C 3 - , 

(Al) 

where, as mentioned in Sec. II, a 2 and a 3 are the second 
(Gauss-Bonnet) and third order Lovelock coefficients. 
The term C\ = R is the Einstein-Hilbert Lagrangian, 
C 2 is the Gauss-Bonnet Lagrangian given by 

C 2 = R lllJy sR lunS - 4+ R 2 , (A2) 

and the third order Lovelock Lagrangian is defined as 

r _ O tdUPCTK TD TlPT I Q D/il^ DCTK, TDPT 

— IHr K aKpr K^ ^ H- OilT ap n VT tir ^ 

_1 OA DFCTK, TD TDP I Q TD TDUVCTK TD 

I Z4xl lt ( j Ki i/ p lL ^ ~T OlL-Ll 

+24 R^™R ail R KV + 16 R^R va R\ 
-URR^R^ + R 3 . (A3) 

In Lovelock gravity only terms with order less than [(n + 
l)/2] (where [x] is the integer part of x ) contribute to 
the field equations, the rest being total derivatives in the 
action. Therefore one should consider n > 6 in the case 
of third order Lovelock gravity so that the effects of third 
order Lovelock gravity appear. 

The (n+l)-dimensional static solution of action (Al) is 
given in the form, which we rewrite for self-completeness, 

dr 2 

ds 2 = - f{r)dt 2 + —- + r 2 dE 2 _ x , (A4) 

where dH^-i is the metric of a (n — l)-dimensional unit 
sphere and the metric function /(r) is given by 

f(r) = 1 — r 2 ip(r). (A5) 

Substituting the electromagnetic potential (2) and the 
metric (A4) into the action (Al) and varying the latter 
with respect to the metric function /(r), one obtains the 
equation of motion as [21] 

a 3 ip 3 (r) + a 2 ip 2 (r) + ip(r) - ^ + r2( ^_ 1} = 0, (A6) 

where, for notational simplicity, we define the following 
coefficients 


a 2 = (n — 3)(n — 2)ct 2 , (A7) 

d 3 = (n ^ 5)(n — 4)(n — 3)(n — 2)03 . (A 8 ) 
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The three solutions of Eq. (A 6 ) are 

+ S+ u5~ 1 , (A9) 

3a 3 

Mr) = --^ + ~(S + uS- 1 )+i^-(S-uS~ 1 ), 

(A10) 

^3( r ) = - i(<5 + M<5 -1 ) - i^y(6 -ud- 1 ), 

(All) 


where 


5 = 


v = 


(v + \] v 2 — w 

9d2<5 3 — 2d2 
545| 


1/3 

5 

1 

2a 3 


u = 


0^2 — 3 3 


j. 2 (n—l) 


9a§ 


m 

rjrtTt 


All of the three roots of Eqs. (A9)-(A11) allow real values 
in appropriate ranges of a 2 and S 3 . For instance, ip\{r) 
is real provided u 3 < v 2 and ip 2 (r) and 1/3 (r) are real for 
u 3 > v 2 . 

In this work, we will only consider A and study some 
properties of the metric function /(r). Note that /(r) = 
1 — r 2 ipi presents an asymptotically flat black hole with 
two horizons provided q < q ex t 7 an extreme black hole 
with one horizon if q = q ext and a naked singularity if q > 
q ex t, where q ex t is related to r ex t through the following 
relationship 


[(« _ 6)53 + („ _ 4 )S 2 r e 2 xt + (n - 2 )r e 4 xt ] - (n = 0 . 

^"ext 

(A12) 

For the special case of = 3a 3 = a 2 , the function 
reduces to 


A(r) = — 


1 - 


1 + 3a 


m 

j.n 


r 2 (ra-l) 


1/3 


(A13) 

In studying wormholes, we consider the region outside 
the horizon r > r^, where is the largest real root of 
/(r) = 0 or = 1 /r 2 . 


Appendix B: Junction conditions in Lovelock gravity 


In order to study these geometries, we first briefly in¬ 
troduce the junction conditions in Lovelock gravity. Let 
us consider two manifolds M+ and M- of M separated 
by a hypersurface E, and denote its two sides by E±. In 
this case, we must add boundary terms to the action so 
that the variation of the action with respect to the metric 
is well-defined. The appropriate boundary terms are [22] 


1 


-^sur — 0 / d *£\/ T 

JdM 


J-2G^>K a 


:(1) T<T ab \ 


K + 2a 2 

+ 3a 3 (P - 2G i * b ) K ab - UR ab J ab + 2 RJ 
-4:R abcd (2 K ac K b K ed - KK ac K bd )\ 1 , (Bl) 
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where J and P are the traces of 

J ab = \(2KK ac K£+K cd K cd K ab -2K ac K cd K db -K 1 2 K ab ) : 

(B2) 

and 

P ab = \{[K A - QK 2 K cd K cd + 8KK cd K d K ec 
5 

-6K cd K de K ef K fc + 3 {K cd K cd f]K ab 
-(4 K 3 4 5 6 7 8 9 10 - 12 KK ed K ed + 8K de K e f Kf d )K ac K£ 
—24:KK ac K cd K de K b + 24K ac K cd K de K ef K bf 
+ (12K 2 - 12K ef K e f) K ac K cd K db } , (B3) 


respectively, where K ab is the extrinsic curvature. 

The tensors G^ b and are respectively the Ein¬ 
stein and Gauss-Bonnet tensors associated with the in¬ 
duced metric 7 ab . Then, the variation of the total action 
/ + Ism with respect to 7 ab provides the surface energy- 
momentum tensor T^. The discontinuous change of the 
surface energy-momentum tensor on the boundary E or 
junction condition can be written as 

Sab = {T+)ab ~ {T-)ab, (B4) 

where the symmetric tensor Tab is given by 


J 


Tab = ^{ R ab - Klab + 2a 2 [3 J ab - J lab - 2G (a c K b)c + 2R ab K - K ab R + 2K cd G cd lab - 2K cd R acbd ] 

+3a 3 {5 P ab - Plab + 2KG ( T + C 2 (K ab + K lab ) + 4 JR ab - 24 J {a c R b)c + 8 K cd R ac R bd - 8 KK a c K b d R cd 
+8KK ab K cd R cd - 8K cd R ab R cd + 16 K {a c R b) d R cd + 16 K (a c K b) d K d e R ce - 8 K ab K c e K cd R de + 6 J ab R 
-8 K (a c R b)c R - 12 J cd R acbd - 4 K cd RR acbd - 8K a c K bc K de R de + 1 6K {a c K c d K e f R b)edf + 16 K {a c K/K de R b)ecf 
+16 K (a c R de R b)dce - 16KK (a c K de R b)dce - 16 K cd R (a e R b)cde - 4 K cd R ac e IR bdef + 8 K cd R c e R aebd 
+8 K cd R c e R adbe - 8K cd R a e jR bedf - 8 K {a c R b) de IR cdef + 8 K a c K b d K e fR cedf - 4 K ab K cd K e IR cedf 
+8 K cd R a e b lR cedf + 2 lab {G { ^K cd + 6 J cd R cd - JR+ 2KK cd K e IR cedf - 4 K c e K cd K fh R dfeh )}. (B5) 
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